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Abstract 

Ballistic deposition (BD) is considered to be a paradigmatic discrete growth model that repre- 
sents the Kardar-Parisi-Zhang (KPZ) universality class. In this paper we question this connection 
by rigorously deriving a formal continuum equation from the BD microscopic rules, which deviates 
from the KPZ equation. In one dimension these deviations are not important in the presence of 
noise, but for higher dimensions or when considering deterministic evolution they are very relevant. 
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I. INTRODUCTION 



Kinetic roughening of nonequilibrium surface growth has been of great interest in re- 
cent years. The kinetic growth processes have been intensively studied via various discrete 
models and continuum equations and exhibit nontrivial scaling behavior 
The surface width W, which is the standard deviation of the surface height, scales as 
W {L,t) ~ L'^g {t/L^), where the scaling function g {u) is constant for u ^ 1 and behaves 
like for m ^ 1. Since the growth is independent of the system size L at the beginning of 
the process, the exponents must obey the scaling relation f3 = a/z. The scaling behavior of 
the growth is characterized by the roughness exponent a, the growth exponent (3, and the 
dynamical exponent z, and these exponents determine the universality class. 

Given a discrete model, it can be very difficult to assign to it a continuum equation that 
describes its dynamics, and vice versa. Nowadays, this kind of problems is usually tackled 
using extensive simulations that determine the critical exponents, by which the universality 
class is determined. A continuum system characterized by the same exponents is then said 
to belong to the same universality class. 

n 

For example, the ballistic deposition (BD) model ^2], which is a discrete model, is believed 



to be described by the famous Kardar-Parisi-Zhang (KPZ) equation [8|. To be more specific, 
we will briefiy present these two models, since they will interest us in this paper. The ballistic 
deposition (BD) model in 1 + 1 dimension (actually on a two-dimensional square lattice) can 
be described as follows. At time t, the height of the interface of site i is hi (t) . We choose 
a random position above the surface and allow a particle to fall vertically toward it. The 
particle sticks to the first site along the trajectory that has an occupied nearest neighbor. If 
no such neighbor exists, it lands on the surface below. Actually, the version just described 
is called NN (nearest neighbor) BD. In another version of the model, also known as NNN 
(next nearest neighbor) BD P], which will interest us in this paper, the particle is allowed 
to stick to a diagonal neighbor as well, as shown in Fig. (0). At time t + 1, a column i is 
chosen at random, and the height hi{t + 1) is then given by 

hi (t + 1) = max {hi_i (t) , h, (t) , h,+i (t)} + 1. (1) 

Not surprisingly, it turns out that both NN and NNN models are well describes by the 
critical exponents of the KPZ equation, as they share the same basic growth mechanism. 
Actually, even more sophisticated versions of the model, including deposition of n particles 



at a time, and a similar model that describes formation of foam (see ref. j6|) show the same 
universal behavior. 
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FIG. 1: Schematic representation of the next-nearest-neighbor bahistic deposition (NNN BD) 
model. A particle falls vertically and sticks to the first site along its trajectory that has an 
occupied nearest neighbor. The particle is allowed to stick to a diagonal neighbor as well. 

As mentioned above, the continuum equation that is believed to capture the essential 
dynamics of the BD models is the famous KPZ equation, given by 

f)h A 

— (f, t) = vV'h (f, t) + - [Whf + (f, t) , (2) 

where h (r, t) is height of the interface above f at time t, and rj (r*, t) is a noise term such 
that 

(r/(f,t)) = 

{'n (f, t) f] (r^, t)) = 2Do6 (f - r^) 6 {t - t') . (3) 

Many proposals for establishing more direct associations between discrete models and 
continuum equations of motion are now present. First, there are phenomenological Ql and 
symmetry arguments that can be very illuminating, but have limited power since 

they cannot connect microscopic quantities to macroscopic ones. Sometimes it is possible 
to map the growth problem onto other models Q| , for which reliable results are known. 
A different approach, based on real-space renormalization-group methods is designed 
to identify the relevant macroscopic constants (for example, diffusion coefficient) from the 
numerical data, once the universality class of the model is already known. 

However, in the last few years most of the efforts in establishing direct connections be- 
tween discrete models and continuum equations have been directed towards developing a 



eneral procedure via formal expansions of discrete equations of motion 



m 0, m Q, 



2G 
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22]. Usually, the derivation of the continuum equation is based 



on regularizing and coarse graining discrete Langevin equations that are obtained from a 
Kramers-Moyal expansion of the master equation. In simple words, transition probabili- 
ties are calculated from the microscopic rules of the model for any given discrete height 
configuration {hi}. These expressions usually contain discrete 9 (Heaviside) and 6 (Dirac) 
functions. But, since the transition probabilities are supposed to be continuous functions 
(so that the expansion that is used there is meaningful) some coarse-graining procedure is 
needed. More specifically, this involves expansions of the form 



9{x) = l + ^Akx' 

k=l 



(4) 



Sometimes a less restricted form 9 (x) 



J2 is used. 

fc=0 



as originally suggested in [l^ 
Other suggestions include 

^(a;) = [1 + tanh(Ca;)]/2, (5) 

where C is an arbitrary positive parameter, with the exact 9 [x) function being obtained 
in the limit C ^ oo |l7[. C is then used in the expansion as an uncontrolled parameter. 
Others ^| use the shifted form 9 (x) = Yim [1 + tanh C {x + a;})]/2 with a G , or a 

modified version using arctan (Cx) 13| or erf [Cx] 20 1 instead of tanh (Cx) (each version 
has its advantages) in eq. 

In some cases the master equation approach is problematic, like in the case of deriving 
the Kardar-Parisi-Zhang (KPZ) equation from the discrete model known as Ballistic 
Deposition (BD) Thus a similar approach, yet more appropriate for that specific case 
was developed. The method is based on dealing directly with the discrete Langevin equation 
rather than with its associated master equation. Still, expansions like 



9{x 



n=l 



(6) 



y=x 



were used. A different, yet closely related Langevin based, approach used the following 
representation of the max function [2^ I21I 1 



max {A, B, C} = lim e In {e^/" + e^^' + e^/"} 



(7) 



in order to go from the discrete BD model to the continuum KPZ equation. Lately 



influenced by the last representation, the Edwards- Wilkinson pj equation was derived from 



a discrete model using 

"g{x+a)/e _|_ 



6 ix) = max {x + a, 0| — max {x} = lim < — In 



(8) 



where a is any constant in the interval (0, 1]. 

In spite of these many new and interesting derivations, one can easily point out three 
main drawbacks of this last approach. First, in many cases the derivation is performed in 
one dimension, where higher dimensions are not discussed at all, or specifically known to 
cause fatal difficulties (see ref. |2l|] for example). 

Second, obviously the mathematics used in the derivation is not so "kosher" . For example, 
an expansion like the one given in eq. Q is problematic because the Heaviside function 
is certainly not analytic around zero. Another example is when taking an expression like 
eq. ((HI) and expanding it for small C - while the limiting procedure that is needed for the 
equality to hold requires C ^ oo. 

Third, since artificial parameters like C and e that cannot always be removed later enter 
the discussion it is not possible to infer the macroscopic quantities (such as the diffusion 
coefficient) from the microscopic rules. 

In this paper, we focus on the ballistic deposition (BD) model, and its believed relation 

nn 

with the Kardar-Parisi- Zhang (KPZ) equation, since it is more difficult to establish [la, |2J| 
than the relation of KPZ to other discrete models such as solid-on-solid (SOS). In what 
follows we will question this relation, and claim that it is not a coincidence that such a formal 
(and "kosher") derivation was not found. In section II we will show that strictly speaking the 
continuum model that describes BD in one dimension is not the KPZ equation, but rather an 
equation with | V/i| instead of the quadratic (V/i)^ term. Then, we will claim that in the case 
of stochastically driven growth, the difference between the two equations is not dramatic, and 
just modifies the value of the coupling constant in the KPZ equation. In section III we show 
(using [3; El) that when the noise is shut down and the dynamics becomes deterministic 
the difference between the BD and KPZ is very important. In section IV, we make a step 
forward towards a derivation of a continuum equation that emanates from BD in higher 
dimensions. However, the equation we find reflects the underlying structure of the discrete 
lattice, in the sense that the equation depends on the directions of the coordinate system 
induced by the lattice. Therefore, the equation we obtain cannot be conclusively related to 
the KPZ equation. At the end, in section V, a brief summary of the results obtained in this 



paper is presented. 



II. DERIVATION OF A CONTINUUM EQUATION IN ONE-DIMENSION 

We begin with the following one-dimensional discrete model 

hi (t + 1) = max {hi^i (t) , hi (t) , hi+i (t)} + r]i (t) , (9) 

where hi (t) is the height of a surface at lattice position i and time t (the time is discrete as 
well). In addition, rji (t) is a white noise term satisfying 

{r],{t)r],{t')) = 2Do6,,,,6t,t'. (10) 
This discrete model is a formal expression for an n-particle NNN (next nearest neighbor) 



BD model (see eq. (21) in ref. 2l|). It is well known that the fact that here we deposit more 
than one particle per unit time is not important for the universal behavior of the model. 
Actually the same discrete model is used to describe formation of foam j^. 
First we represent the max{} operator using Heaviside functions 

max {a, b} = a6 {a — b) + bO {b — a) . (11) 

Thus, for three arguments we also get 

max {a, b, c} = aO {a - b) 6 {a - c) + bO {b-a)e{b-c) + c9 {c-a)9{c-b). (12) 

Now we represent the Heaviside function using the sign function 

^(a;) = ^ + ^sgn(x). (13) 

Applying these representations to eq. (jHI) we get 

max {hi^i (t) , hi (t) , hi+i (t)} 
= (t) i [1 + sgn (/ii_i (t) - hi (t))] [1 + sgn (/i,_i (t) - (t))] 
+ hi (t) i [1 + sgn {hi (t) - hi^i (t))] [1 + sgn {hi (t) - hi+i (t))] 
+ hi+^ it) \ [1 + sgn (/li+i it) - (t))] [1 + sgn {hi+i it) - hi it))] . (14) 
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Thus, after some simple algebra 



hi (t + 1) - hi {t) = I [hi^i (t) - 2h (t) + hi+i (t)] - \hi (t) 

+ \ [hi it) - it)] sgn {hi (t) - (t)) 

+ i [hi+i (t) - hi^i (t)] sgn {hi+i (t) - h^^l (t)) 

+ i [hi+i (t) - hi (t)] sgn (t) - hi (t)) 

+ \hi_i (t) sgn {hi (t) - hi_i {t)) sgn (/ij+i (t) - hi^i {t)) 

- \hi (t) sgn {hi (t) - hi^i (t)) sgn (/i^+i (t) - /i^ (t)) 
+ (t) sgn (t) - (t)) sgn {hi+i (t) - (t)) + 7]i (t) . (15) 

The last expression is written in such a way that identifying discrete derivatives is easy. 
Therefore, denoting the spatial and temporal increments by Ax and At, and using the 
simplifying fact sgn (ax) = sgn (x) (when a > 0) we get 

^2 r /FthM r)27, r)h 1 

(16) 



dh (Ax) 



1 + sgn^ — 



dx I 



d^h , , Ax 
At 



dx 



Note that sgn^ (x) 



X ^ 



. Thus, for a strictly non smooth surface (i.e. a surface 



1 

X = 

that is not fiat), almost everywhere we can use the replacement sgn^ {^^^ ~ ^ order 
to further simplify eq. (fTH|) . In addition, since in the discrete model we actually took 
Ax = At = 1 we get the final result 

dh 



dh , , 1 d'^h . , 



dx 



+ r/(x,t). 



(17) 



This equation is a result of a straightforward exact coarse-graining of eq. ©. Thus, one 
can be convinced that BD (at least this n-particle version of it) does not lead exactly to 
the KPZ equation. This fact might explain the difficulties encountered in the past while 
deriving the KPZ equation from BD. 

At this point, it is interesting to wonder what is the connection between BD and the KPZ 
equation? After all, for the last two decades BD is considered to be a paradigmatic discrete 
growth model that belongs to the KPZ universality class (see for example the approach taken 
by Barabasi and Stanley p to present the field of surface growth). In addition, extensive 



simulations 



26 



27 M 



30| aimed at extracting the KPZ critical exponents used 



BD, and found critical exponents that are consistent with other known results. However, one 
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that 



important point should be emphasized. Most of the simulations |25 . 
used BD were done in one and two dimensions, while only one of them (ref. |2^) studied BD 
in three dimensions. This point is important because in one and two dimensions only a rough 
phase is possible, while for higher dimensions a phase transition between a smooth phase and 
a rough phase becomes possible. The last remark will be a starting point for the discussion 
in section IV, where the behavior of BD in higher dimensions in discussed. In the following, 
we will discuss, however, the equivalence of the one dimensional continuous model (eq. IT7|) 
we derived from the BD deposition discrete model. A support to that equivalence is already 
found in ref. j3| that studied stochastically driven Langevin equations (of the KPZ type) 
with a generalized |V/?.|^ nonlinearity. They concluded that no-matter what value of fi was 
taken, KPZ behavior was revealed, as long as the noise term is present (the case of no noise 
is interesting by itself, and will be discussed in section III below) . This result should not be 
a surprise nowadays as it is consistent with the well known symmetry argument [3, 0, 0, Q| 
that is often used to justify the KPZ equation -namely since the generalized nonlinearity 
|V/?.|^ breaks the h —>■ —h symmetry it should belong to the KPZ universality class. Thus, 
this fundamental observation of Amar and Family [2^ establishes the link between eq. fll7|l 
and the KPZ equation, by simply identifying /i = 1 in that equation. Hence, the route 
between BD and KPZ in one dimension is now understood, and the common knowledge 
about BD is justified. 

Now, it would be interesting to complete the picture and to extract the macroscopic 
coefficients that describe the continuum KPZ equation that is related to the BD model. 
First, the diffusion coefficient can be inferred easily from equation (fTTj) . Then, the noise 
amplitude is unaffected by the coarse graining. Finally, we would like to identify the coupling 
constant A in the KPZ equation (j21). However, since the nonlinear term in eq. ()17p is not 
in a KPZ form, A cannot be just read from the equation. Thus, we push further the idea 
presented in the previous paragraph, and make the following identification 

\Vh\ \Vhf 



(IV/^I) (\Vh 



|2\ ' 



(18) 

where (■ ■ ■) means steady-state averaging. The last equation leads to|V/;-| ^ (^\vhf) 
Plugging this estimate into eq. (fTTj) . gives the prediction A/2 ~ ^jl^^^^ for the coupling 
constant of the KPZ equation Q. The meaning of this expression is that the KPZ coupling 
constant that describes the surface that grows under BD is obtained by calculating the 
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two quantities (|V/?.|) and (j'Vhf'^ in a BD simulation (in steady state) and plugging the 
obtained averages into the expression given above for A/2. 

The last prediction has been checked numerically by us. We grew a one dimensional 
interface on a discrete lattice with linear size L = 1024 using the microscopic rules of BD 
(given by eq. Q above). Using this simulation we found (|V/;,|) = 0.4705 and <^|V/i|^^ = 
0.4639 that leads to the prediction A ^ 2.17. Using an inverse method - that is "guessing" 
the best macroscopic parameters that recover the same growth process (see refs. ^ ^J), 
we found A = 2.27, that supports the rough estimate of A we suggested above. 

Up to here we focused on the common characteristics of BD and KPZ. However, as will be 
seen in the following two sections, if the noise is turned off, or if we go to higher dimensions, 
significant differences between the two models can be found. 

III. DETERMINISTIC FLATTENING 

In this section we discuss the dynamics of the discrete BD model and the continuum 
model with |V/i|^ non-linearity for ^ = 2 (KPZ) and for = 1, when the noise term that 
appears in them is eliminated. This limit mimics the physical scenario when the deposited 
materials stop falling, and the interface relaxes to a flat surface, due to its diffusive term. We 
think that it is worth studying such physical situations because it emphasizes the difference 
between the dynamics of BD and KPZ already in the one-dimensional case. 

The starting point of this discussion is the observation of Krug and Spohn [2^ regarding 
a deterministic continuum equation of the general form 

dh 

^ = uV'h + g\Vhr, (19) 

(where /i > 1) that describes the smoothing of an initially rough surface under deterministic 
growth. They argued that if the initial surface has roughness exponent a, the scaling relation 

Zd = min {2, /x (1 — a) + a} (20) 

should hold, where Zd is a dynamic exponent for deterministic evolution, which similar to 
the stochastically driven case, determines the time scale of the decaying width (see Fig. |2I). 

n 

Influenced by this work, Amar and Family :24] made a systematic study of such surface 
growth equations with a generalized nonlinearity and found good agreement between the 
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scaling relation and the numerical integration of equation ()19|). We also performed a nu- 
merical integration, using the same integration scheme as in ref. Q] (the coefficients v and 
g were taken from eq. ()17|1 and using the result of the previous section) for the cases yU = 1 
(eq. (fTTjl ) and /i = 2 (KPZ) in 1+1 dimension. We also used a rough surface with a = | as 
an initial condition. A scaling plot of our results for the surface width W {L, t) for fi = 1 
versus t/L^ with Zd = 1 is given in Fig. Efa), and a similar scaling plot for /i = 2 with = ^ 
is shown in Fig. I2tb). 




FIG. 2: Scaling plot for the deterministic equation (|19jl for L = 256,512 and 1024 with (a) = 1 
using Zd = I and with (b) /i = 2 using = |. 



As one can see, the scaling is good in both cases, in agreement with the scaling relation 
fl2()|l and therefore consistent with refs. 3, • Thus, the difference between the deter- 
ministic dynamics of the KPZ equation and that of eq. p7|l is now self evident, and not 
new. Two questions should be asked here. Our derivation of eq. ()17|) does not depend on 
the presence of noise. Therefore, we might expect that even when deterministic dynamics 
is considered discrete BD will be equivalent to the fi = 1 case rather than to the fi = 2 case 
(KPZ). The first question is whether this is indeed the case. The second question regards 
the predicted values for the diffusion coefficient (namely z/ = |) and the coupling (i.e., the 
coefficient of the | V/i]'^ term - namely g = 1) for the n-particle NNN BD model that we read 
from equation ()17p. The question is whether this prediction is a sensible one? 

In order to answer these questions we implemented the NNN BD model (whose sticking 
rules are given by eq. Q) on a two dimensional square lattice (1 + 1 dimensions), and 
measured the surface width W {L, t) of the smoothing surface, using a rough surface with 
a = ^ as an initial condition (just like the continuum models above). We present the results 
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in Fig. El 




t/L 



FIG. 3: Scaling plot for the deterministic NNN BD model given by eq. © for L = 256, 512 and 
1024 using = I. 

As can be seen in Fig. El the deterministic dynamics of the discrete model is evidently 
the one we have seen for the = 1 continuum equation, since the data collapse in the scaling 
plot with 2;^ = 1 is good. Moreover, Fig. Efa) and Fig. 01 look pretty much the same. The 
main difference between the two is the rounded tail in Fig. Efa) around t/L = 0.5 compared 
to the sharp transition in Fig. El This difference can be accounted for by the fact that the 
numerical integration of the continuum equation used a small time increment (At = 0.05) 
while the discrete growth model used a much larger one (At = 1). This implies a better 
temporal resolution around t/L = 0.5 by the numerical integration when compared to the 
simulation. I addition, the similarity between the two scaling plots means that the predicted 
macroscopic quantities are consistent with the results obtained from the simulation. 

In this section we showed that a difference between BD model and KPZ equation can be 
seen, when looking at their deterministic dynamics. We also showed that eq. ()17|) captures 
the dynamical behavior of BD even in the deterministic regime. 

IV. RESULTS FOR HIGHER DIMENSIONS 

In this section we generalize the formal derivation given in section II to higher dimensions. 
As will be seen, the derivation is not exactly the same, and the resulting continuum equation 
is not a simple generalization of the one dimensional equation ()17|1 . 
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We begin with the n-particle NNN BD model in d dimensions (discretized on a cubic 
hyperlattice) 



h,^ it + l) = max {h,^ {t) , (t) , (t) , . . . , K^+x^ it) , h,^-^^ (t)} + 77^- (t) , (21) 

where Xi IS cL unit vector in the i'th direction. 

Like in section II, using equations (fT^ - (fT!?jl . coarse graining (space) and using the sim- 
plifications sgn [ax] = sgn (x) (for a 7^ 0) and sgn^ (^) ~ -'- (which is again true almost 
anywhere) we get 

2'%^{t + 1) = 2X: h,^,^ it) [1 + sgn {§-)] n [1 + sgn (g - §-)] [l + sgn {§- + §- 



i=l 



+ 



+h,.,^ it) [1 - sgn {§-)] n [1 - sgn {§- - §-)] [1 - sgn {§- + §- 
Using 1 ± sgn (x) = 26 (±a;) we get 

A, (i + 1) = i h,,,. (t) (ft) n [» ( ft - ft) « (ft + i 



+ *(«)■ (22) 



1=1 



+ 



(«) " (-ft) n K-fe +£)"(-£- 1 



+ r/,-.(t), (23) 



where r)^?(t) = 2 '^'^r]f^{t). Notice that 

e{a)e{a + x)e{a-x) = e (a) ^ (a - |x|) = 6 (a) 6 (|a| - |x|) 



(24) 



thus 



d d 



i=i j=i 



dh 



dxi 



dh 



dx-j 



hr+^. it) ^ (fe) + hr-^. it) ^ {-§-)] + Vrit) . (25) 



Using once again the relation 6* (x) = | [1 + sgn (x)], and reorganizing the last equation lead 
us to 

d 



hrit+i) = iEn^( 



i=l j=l 



dh 




dh 




dxi 




dXj 


) 



X 



K-+^^ (t) - 2hr- {t) + K-_^^ (t) + 2gsgn [§-) + 2h^ (t)] + fj,. it) , (26) 
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Coarse graining in both space and time gives us the final result 



dh 



i=l 



2dx] 



+ 



dh 



dxi 



n^( 



dh 



dxi 



dh 



) +v{f,t). 



(27) 



This equation might seem messy for a second. However, it has a very simple structure. 
First, for d = 1 it collapses trivially to eq. ()17|) . Then, in higher dimensions, the expression 



n H 

1=1 ^ 



3 



dh 




dh 




dxi 




dXj 


) 



picks the direction, m, in space along which the gradient is maximal, and 



the considered part of the local growth rate has a contribution from that direction alone, 

1 d^h 

2 dxi 



+ 



dh 



dXr, 



Actually, taking a glance at the original discrete model, this result is not such 
a surprise since both formulations contain this "maximal" growth ingredient in them. 

Actually, the equation we obtain still refiects the directions that are induced by the 
underlying lattice. Therefore, strictly speaking, this is not a proper continuum description. 
Moreover, it is not clear whether such a proper continuum description exists, and how can 
it be derived. Therefore, the relation of this equation to the d-dimensional KPZ equation 
for d > 2 is also not evident. We leave these questions for future research. 



V. SUMMARY AND CONCLUSIONS 



In this paper we discussed the connection between a discrete growth model, the next- 
nearest-neighbor ballistic deposition, and a continuum equation, the Kardar Parisi-Zhang 
equation. It has been believed for the last two decades [^, that BD is a paradigmatic 
discrete model for the KPZ universality class; however a formal derivation was lacking. 
In this work, we show that the absence of a formal derivation is not accidental, but rather 
refiects significant differences between the continuum equation that describes the BD model, 
and the KPZ equation. This difference turns out to be mild in one-dimension in the presence 
of noise, but crucial when discussing the deterministic dynamics. The relation between BD 
and KPZ is also questioned when going to higher dimensions. 

The main advantage of our approach is that it is direct and allows to estimate of the 
macroscopic quantities of interest (such as the diffusion coefficient) from the microscopic 
rules. 

The main disadvantage of our approach is that it does not give a general prescription for 
deriving continuum equations from discrete models. However, we do show that the option 
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of giving a rigorous formal derivation should not be overlooked. As we show in this paper, 
the usage of unjustified expansions (such as those presented in eqs. (jH)-® above) might 
sometimes lead to wrong conclusions. 

Two open questions are left for future work. First, whether the BD model in d dimensions 
has a proper continuum description, which does not depend on the discrete lattice on which 
it is defined, and how can one derive this equation. Second, whether the BD model belongs 
to the KPZ universality class in dimensions higher than one. 
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